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Abstract 

We introduce two basic invariant forms which define generic surface in 
3-space uniquely up to Lie sphere equivalence. Two particularly interesting 
classes of surfaces associated with these invariants are considered, namely, 
the Lie-minimal surfaces and the diagonally-cyclidic surfaces. For diagonally- 
cyclidic surfaces we derive the stationary modified Veselov-Novikov equation, 
whose role in the theory of these surfaces is similar to that of Calapso's 
equation in the theory of isothermic surfaces. Since Calapso's equation it- 
self turns out to be related to the stationary Davey-Stewartson equation, 
these results shed some new light on differential geometry of the stationary 
Davey-Stewartson hierarchy. Diagonally-cyclidic surfaces are the natural Lie 
sphere analogs of the isothermally-asymptotic surfaces in projective differen- 
tial geometry for which we also derive the stationary modified Veselov-Novikov 
equation with the different real reduction. 

Parallels between invariants of surfaces in Lie sphere geometry and recip- 
rocal invariants of hydrodynamic type systems are drawn in the conclusion. 
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1 Introduction 



Lie sphere geometry dates back to the dissertation of Lie in 1872 After that 
the subject was extensively developed by Blaschke and his coworkers and resulted in 
publication in 1929 of Blaschke's Vorlesungen iiber Differentialgeometrie ||, entirely 
devoted to the Lie sphere geometry of curves and surfaces. The modern multidi- 
mensional period of the theory was initiated by Pinkall's classification of Dupin 
hypersurfaces in E 4 ||, jp. We refer also to Cecil's book || with the review of the 
last results in this direction. Since most of the recent research in Lie sphere geometry 
is concentrated around Dupin hypersurfaces and Dupin submanifolds, the general 
theory of Lie-geometric hypersurfaces seems not to be constructed so far. The aim 
of this paper is to shed some new light on Lie sphere geometry of (hyper)surfaces 
and to reveal its remarkable applications in the modern theory of integrable systems. 

Let M 2 e E 3 be a surface in the 3-dimensional space E 3 parametrized by the 
coordinates R 1 , R 2 of the lines of curvature. Let k 1 , k 2 and gudR l2 + g 22 dR 22 be the 
principal curvatures and the metric of M 2 , respectively. In sect. 2 we introduce the 
basic Lie sphere invariants of the surface M 2 , namely, the symmetric 2- form 



WW 



dR l dR 2 (1.1) 



(A; 1 — k 2 ) 2 

and the conformal class of the cubic form 

WgndR 13 + d 2 k 2 g 22 dR 23 J (1.2) 

di = t^-, which define "generic" surface M 2 uniquely up to Lie sphere equivalence. 
We recall that the group of Lie sphere transformations in E n+1 is a contact group, 
generated by conformal transformations and normal shifts, translating each point 
of the surface to a fixed distance a = const along the normal direction. Conformal 
transformations and normal shifts generate in E n+l a finite-dimensional Lie group 
isomorphic to SO(n + 2, 2). Lie sphere transformations can be equivalently charac- 
terized as the contact transformations, mapping spheres into spheres and preserving 
their oriented contact. In the implicit form objects (1.1) and (1.2) were contained 
already in Quadratic form (1.1) gives rise to the Lie-invariant functional 

dtfdR 2 (1.3) 



(A; 1 - k 2 ) 2 

the extremals of which are known as minimal surfaces in Lie sphere geometry (||, 
§ 94). Other possible representations of functional (1.3) and its relationship to 
the integrable hierarchy of Davey-Stewartson type are discussed in sect. 3. In 
sect. 4 we investigate the so-called diagonally-cyclidic surfaces ("diagonalzyklidische 
flachen" in the terminology of [[|, p. 406), which can be characterized as surfaces M 2 
possessing parametrization R 1 , R 2 by the coordinates of lines of curvature such that 
the cubic form (1.2) becomes proportional to dR lS + dR 2 . This class of surfaces is 
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a straightforward generalization of isothermic surfaces in conformal geometry. It is 
demonstrated that for diagonally-cyclidic surfaces the Lie sphere density U defined 
by 

2 _ d\k x d 2 k 2 

satisfies the stationary modified Veselov-Novikov (mVN) equation 

dfU + 3VdiU + \UdxV = dlU + 3Wd 2 U + \Ud 2 W 
d x W = d 2 {U 2 ) 
d 2 V = d 1 {U 2 ) 

which in the theory of diagonally-cyclidic surfaces plays a role similar to that of the 
Calapso equation [|| 

2 / did 2 u\ 2 ( did 2 u\ 1 2 

in the theory of isothermic surfaces. Since Calapso's equation itself turns out to 
be related to the stationary Davey-Stewartson (DS) equation (see sect. 4), this pro- 
vides remarkable differential-geometric interpretation of the stationary flows of DS 



hierarchy and gives new insight into the recent results of || - |15j relating DS hi- 
erarchy with conformal geometry. The details of derivation of mVN equation for 
diagonally-cyclidic surfaces are included in the Appendix. 

It is quite remarkable that in projective differential geometry there also exists 
a class of surfaces (the so-called isothermally-asymptotic surfaces, or ^-surfaces) 
governed by the stationary mVN equation 

PXXX ~ 3Vp x - \pV X = Pyyy ~ 3W Py ~ \pWy 



W x = (p 2 

Vy = (P 2 



which can be reduced to that presented above by a complex change of variables 
p — > ill, W — ► —W, V — > —V. These are surfaces, for which the Darboux cubic 
form is proportional to the sum of pure cubes dx 3 +dy 3 in the appropriate asymptotic 
parametrization x, y (particular examples are affine spheres related to the Tzitzeica 
equation). This observation reflects the duality between projective and Lie sphere 
geometries due to the Lie's famous line-sphere correspondence (see sect. 4). 

In sect. 5 we introduce Lie sphere invariants of higher dimensional hypersurfaces 
M n G E n+1 , namely the symmetric 2- form 

h.ih.3 

and the conformal class of the cubic form 

£%iio/ 3 (1.5) 
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defining "generic" hypersurface uniquely up to Lie sphere equivalence. Here k l 
are principal curvatures, uj % are principal covectors, Yj gu^ 1 is the first fundamental 
form and coefficients k\ are defined by the expansions dk % = k^uji (we emphasize that 
hypersurface M n of dimension n > 3 does not necessarily possess parametrization by 
the coordinates of lines of curvature). Objects (1.4) and (1.5) are the Lie-geometric 
analogs of the second fundamental form and the Darboux cubic form in projective 
differential geometry of hypersurfaces. 

In sect. 6-9 interrelations between Lie sphere invariants and reciprocal invariants 
of hydrodynamic type systems 

u\ = v){u)u> x i i,j = l,...,n (1.6) 

are discussed. We recall that reciprocal transformations are transformations from 
x, t to the new independent variables X, T defined by the formulae 

dX = B(u)dx + A(u)dt 
dT = N(u)dx + M{u)dt 

where Bdx + Adt and Ndx + Mdt are two integrals of system (1.6). Reciprocal 
transformations originate from gas dynamics and were extensively investigated in 
[ |l~6|j , [|i~7f . In |nj, ||19[ we introduced reciprocal invariants, defining a hydrodynamic 



type system uniquely up to reciprocal equivalence. The summary of these results in 
the 2-component case is given in sect. 6. In sect. 7-8 we recall the necessary infor- 
mation about Hamiltonian systems of hydrodynamic type and describe the general 
construction of [20|, [21], relating Hamiltonian systems (1.6) and hypersurfaces in 



E n+l . The main property of this correspondence is its " equivariance" in the sense 
that Lie sphere transformations of hypersurfaces correspond to "canonical" recip- 
rocal transformations, that is, to those reciprocal transformations which preserve 
the Hamiltonian structure. In this approach Lie sphere invariants of hypersurfaces 
correspond to reciprocal invariants of hydrodynamic type systems, providing thus 
their differential-geometric interpretation. 

In sect. 9 we write down reciprocal invariants of n-component systems for arbi- 
trary n > 3 since they differ from those in case n — 2. 



2 Invariants of surfaces in Lie sphere geometry 

In UJ], p. 392 Blaschke introduced the Lie-invariant 1-forms uj 1 , uj 2 (dtp, dijj in 
Blaschke's notation) which assume the following form in the coordinates R 1 ,^ of 
the lines of curvature: 



(2.1) 
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Remark 1. In order to check Lie-sphere invariance of the 1-forms uj 1 ,uj 2 it is 
sufficient to check their invariance under the inversions and normal shifts, which 
can be verified by a direct calculation. Moreover, forms (2.1) do not change if the 
principal curvatures k % are replaced by the radii of principal curvatures w i = h and 
the first fundamental form gn by the third fundamental form Gu = k l g^. 

Remark 2. Similar invariant 1-forms arise in the Mobius (conformal) and La- 
guerre geometries, which are subcases of the Lie sphere geometry. In the Mobius 
geometry we have the invariant 1-forms 

^ :dR\ J^dR 2 



k 1 — k 2 ' k 2 - k 1 
and the invariant quadratic form 

(k 1 -k 2 ) 2 (g ll dR l2 + g 22 dR 22 ), 

while in the Laguerre geometry they are 

9lWl dR\ 4^-rdR 2 



w 1 — w 2 w 2 — w 1 



and 

(w 1 - w 2 ) 2 {GudR l2 + G 22 dR 22 ), 

respectively (see @, |28|, @). 

As far as uj 1 and u 2 are invariant under Lie sphere transformations, so do the 
quadratic form 



-uW = dlkl9 ^l dtfdR 2 



(k 1 - k 2 ) 



2 



and the qubic form 



^ - co 23 = m dl ^ 2k l (d 1 k 1 g 11 dR 13 + d 2 k 2 g 22 dR 2 % (2.2) 

[k — K ) y/giig 22 

giving rise to (1.1) and (1.2), respectively. The reason for introducing these objects 
is their additional invariance under the interchange of indices 1 and 2 so that they 
play a role similar to that of "symmetric functions of the roots of polynomial" in 
the Viete theorem. Hence they define tensors which can be effectively computed in 
an arbitrary coordinate system (without solving algebraic equations). 

In 0, §85 it is proved that up to certain exeptional cases a generic surface in 
E 3 is determined by the corresponding 1-forms uj 1 ,oj 2 uniquely up to Lie sphere 
transformations. Since we can reconstruct uj 1 ,uj 2 from the quadratic form (1.1) and 
the conformal class of the cubic form (1.2) (the multiple in (2.2) is not essential), 
we can formulate the following 
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Theorem 1. A generic surface M 2 e R 3 is defined by the quadratic form 

d k}d k 2 
(Jfci - I 2 ) 2 

and the conformal class of the cubic form 

d^gndR 13 + d 2 k 2 g 22 dR 23 

uniquely up to Lie sphere transformations. 

The vanishing of the qubic form is equivalent to the conditions d\k l = d 2 k 2 = 
which specify the so-called eyelids of Dupin. We recall that the vanishing of Dar- 
boux's cubic form in projective geometry specifies quadrics, which are thus projective 
"duals" of eyelids of Dupin. 

Remark. The principal directions of the surface M 2 can be characterized as 
the zero directions of quadratic form (1.1). On the other hand, they are exactly 
those directions, where cubic form (1.2) reduces to the sum of pure cubes (without 
mixed terms). It should be pointed out, that any cubic form on the plane can 
be reduced to the sum of cubes, and the directions where it assumes the desired 
form are defined uniquely. In projective differential geometry of surfaces these are 
asymptotic directions of the Darboux cubic form. 



3 Minimal surfaces in Lie sphere geometry 

Lie-minimal surfaces are defined as the extremals of Lie- invariant functional (1.3) 

d^dik 2 



(k 1 - k 2 y 

which is a natural analog of the Willmore functional 

/ Jik'-k'f^gT^dR'dR 2 
in conformal geometry and the invariant functional 

J J {w 1 - w 2 ) 2 yjG n G22 dR^R 2 
in the Laguerre geometry. Due to the obvious identity 



(k 1 - k 2 ) 2 (k 1 - k 2 ) 2 {k 1 -k 2 J 

we see that for compact surfaces with k 1 ^ k 2 (for instance, immersed tori) func- 
tional (1.3) coincides with the functional 

/ / w^w dRldR2 =-JJ ab dRldR2 ' (3 - 1} 
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where we introduced the notation a = Jfe^i , b = ■ We recall that in terms 
of the coefficients a and b the Peterson-Codazzi equations of the surface M 2 can be 
written as follows 

d 2 In v /^7 = a, di In yfg^ = b 

while the equation for the radius- vector r assumes the form 

did 2 r = adif+ bd 2 f (3.2) 

manifesting the fact that the net of the lines of curvature is conjugate. 
Introducing the rotation coefficients /3i 2 , (3 2 \ by the formulae 




we can rewrite our functional in the form 

J J fapn dR l dR 2 

which has the meaning of the "integral squared rotation" of the surface M 2 . 

Written in the form (3.1) our functional is closely related to the simplest quadratic 
conservation law of the (2+l)-dimensional hierarchy of the Davey-Stewartson (DS) 
type. To make it clear we recall the construction of which defines the DS-type 
flow on conjugate nets in E 3 . Let M 2 be a surface parametrized by conjugate coordi- 
nates R},R 2 . The radius-vector r of such surface satisfies the equation (3.2) (at the 
moment we do not assume that our conjugate net is the net of lines of curvature). 
Let us define evolution of M 2 with respect to the "time" t by the formula 

ft = ad\ 2 f + (3d2 2 f + pdif + qd 2 f, (3.3) 

where a, (3 = const. The compatibility conditions of (3.2) and (3.3) give rise to the 
integrable system 

at = /3d 2 2 a — ad 2 a + 2adi(ab) + /3<9 2 (a 2 ) + pd\a + d 2 (qa), 
b t = a di 2 b - (3d 2 2 b + 2(3d 2 (ab) + ad^b 2 ) + qd 2 b + d x {pb), 

(3.4) 

d 2 p + 2ad\a = 0, 
d iq + 2(3d 2 b = 0, 

whose relationship to the DS system was clarified in In fact system (3.4) is a 
linear combination of two simpler systems, corresponding to the choices (a — 1, f3 — 
0, q — 0) and (a = 0, (5 = l,p = 0), namely 

a t = —d\ 2 a + 2di(ab) + pdia, 

6 t = d 1 2 6 + 9 1 (6 2 )+9 1 (p6), (3.5) 
d 2 p + 2d x a = 0, 
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and 

at = d 2 2 a + d 2 (a 2 ) + d 2 (qa), 



k = -d 2 2 b + 28 2 {ab) + qd 2 b, (3.6) 
d x q + 2d 2 b = 0, 

respectively. Both systems (3.5) and (3.6) commute according to the general dis- 



cussion m 



and possess the quadratic integral 
j J ab dR l dR 2 

coinciding with (3.1). We emphasize, however, that the nets of lines of curvature are 
not preserved (in general) by these t-evolutions, although always remain conjugate. 



4 Cyclidic curves and diagonally-cyclidic surfaces 

With any surface M 2 we associate a 3-web of curves (that is, three 1-parameter 
families of curves) formed by the lines of curvature and cyclidic curves ("zyklidische 
kurven" in the terminology of Blaschke 0, §86) which are the zero directions of 
cubic form (1.2). In view of formula (2.2) the curves of this 3-web can be defined in 
terms of 1-forms (2.1) by the equations 

£^ = 0, uj 2 = 0, uj 1 -uj 2 = 0, (4.1) 

respectively. Cyclidic curves naturally arise in the attempt to find those eyelids of 
Dupin, which are the "best" tangents to a given surface M 2 at a given point (see 
H, §86 for the details). These curves are the natural analogs of the Darboux curves 
in projective differential geometry. Let us compute the connection 1-form uj of the 
3-web (4.1), that is, the 1-form which is uniquely determined by the equations 



duj 1 = uj A uj 1 , du 2 = uj A uj 



2 



[see 
in 

CO 



27| for the introduction in web geometry). A direct computation results 

H^ + i&SO dR1 + 5 + Mr) dR'+ 

(4.2) 

1/7 In d 1 k 1 d 2 k 2 

Since both uj 1 ,uj 2 are invariant under Lie sphere transformations, so does the con- 
nection 1-form uj. From (4.2) it immediately follows that the curvature form duj of 
3-web (4.1) is given by 

duj = -dQ, 
3 ' 
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where 

( chW , d 1 k 1 \ . f ftfrfci g 2 fc 2 \ , ff2 , . 

As we will see in sect. 6 the object analogous to (4.3) arises in the theory of reciprocal 
invariants of hydrodynamic type systems. 

An interesting class of diagonally-cyclidic surfaces (" diagonalzyklidische flachen" 
in the terminology of @, p. 406) is specified by the requirement, that 3-web (4.1) is 
hexagonal or, equivalently, has zero curvature: 

duj = —dQ = 0. 
3 

In this case there exist coordinates R l ,R 2 along the lines of curvature (note, that 
we have a reparametrization freedom R 1 — > (p l (R 1 )), where uj 1 ,uj 2 assume the form 

u 1 =pdR 1 1 u 2 = -pdR 2 

with nonzero common multiple p. Since in these coordinates the cubic form uo l3 — to 23 
is proportional to dR} 3 + dR 23 , these surfaces are the natural Lie-sphere analogs of 
isothermic surfaces in conformal geometry (we emphasize, that the class of isother- 
mic surfaces is not invariant under the full group of Lie sphere transformations). 

Equations governing diagonally-cyclidic surfaces can be easily written down as 
follows. Since the cubic form 



is proportional to 
we can put 



d^gndR 13 + d 2 k 2 g 22 dR 23 
dR 13 + dR 23 , 



e 2p e 2p 



011 = TTTT, #22 



d X k V * d 2 k 2 ' 

Inserting this representation in the Gauss-Peterson-Codazzi equations, we arrive at 
the following system for k 1 , k 2 , p: 

diP = b+ l 2 9 ^ 

d2P = a + l^, (4.4) 

dik 1 (dxb + |^ In |g) + d 2 k 2 (d 2 a + %d 2 In |g) + fc W = 0, 

where a = b = 

Remark. System (4.4) is a Lie-sphere analog of the system 

dip = b, 

d 2 p = a, (4.5) 
d\p + d 2 p + k 1 k 2 e 2p = 0, 
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describing isothermic surfaces in conformal geometry. The integrability and dis- 
cretization (see ||, 0) of system (4.5) are based on the 50(4, l)-linear problem, 
which comes from the following geometric fact: all isothermic surfaces possess Riba- 
cour transformations preserving the metric up to a conformal factor. Moreover, the 
spectral parameter is due to the following scailing symmetry of system (4.5): 

R 1 = -R\ R 2 = -R 2 , k 1 = ck\ k 2 = ck 2 , c = const, 
c c 

We recall also that system (4.5) can be rewritten as the single fourth-order Calapso 
equation H] 

for the conformal factor u = e p (k 1 — k 2 ). Another possible approach to the integra- 
bility of isothermic surfaces is based on the relationship of the Calapso equation to 
the DS-II equation 

iu t + u xx -u yy + uv = 0, 

Uxx 4" v yy \ u \xx \ U \yy 

which in the stationary case assumes the form 

UXX - Uyy + UV = 0, 

I 1 2 112 

U X X + Uyy = \U\ XX \U\yy 

Excluding v we arrive at the fourth-order equation with respect to u 

A (^^L) + \u\ 2 xx -\u\ly = 



coinciding with (4.6) after the transformation R 1 = x + y, R 2 = x — y and the 
reduction u = u. 

Similar approaches can be applied to system (4.4). Indeed, all diagonally-cyclidic 
surfaces possess Ribacour transformations preserving the cubic form up to a confor- 
mal factor (H], p. 420). Moreover, system (4.4) possesses a similar scaling symmetry 

R} = -R\ R 2 = -R 2 , k 1 = <?k x , k 2 = c 3 k 2 , 
c c 

which is responsible for the spectral parameter. We hope to develop this geometric 
approach elsewhere. Another approach to the integrability of diagonally-cyclidic 
surfaces is based on their remarkable relationship to the mVN equation 

U t = U XXX - Uyyy + 3U X V ~ 3UyW + | UV X ~ \UW y 

W x = (U 2 )y 



Vy = (U 



2 \ 

X 
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introduced in which is the third-order flow in the DS-I hierarchy. In the Ap- 



pendix we demonstrate, that the Lie sphere density U defined by 

2= d 1 k 1 d 2 k 2 
(k 1 - k 2 ) 2 

satisfies in case of diagonally-cyclidic surfaces the stationary mVN equation 

U XXX + 3U X V + \UV X = Uyyy + WyW + \\JW y 
W X = (U 2 )y 

Vy = (U 2 ) x 

(here d\ = d x , d 2 = d y ). Note that U 2 is a conserved density of the mVN equation. 
A passage from equations (4.4) to he mVN equation requires quite complicated 
calculations which were performed with Mathematica (see the Appendix). 

Particular solutions of the stationary mVN equation can be obtained by the 
ansatz 

_2C7^ 1 /C^ 
3 U 3\U 



3 U 3\U 
where U satisfies the Tzitzeica equation 

(\nU) xy = —U 2 + — , c = const. 

In this case the first equation is satisfied identically, while the last two become just 
two conservation laws of the Tzitzeica equation. 

We have demonstrated that the stationary flows of the DS hierarchy have a 
natural interpretation within the contexts of conformal and Lie sphere geometries. 



This observation agrees with the results of 0, |10| , fill , [[L3|], [|Tj], |To| where dy- 
namics of surfaces, induced by the odd flows of DS-II hierarchy, was investigated. 
In particular, it was argued that the integrals of DS-II hierarchy define conformally 
invariant functionals and the corresponding stationary flows define certain confor- 
mally invariant classes of surfaces. We hope that results presented above contribute 
to these investigations. We emphasize also, that stationary points of some of the 
DS-integrals are probably invariant under the full group of Lie sphere transforma- 
tions, rather then just conformal group (for example, the class of diagonally-cyclidic 
surfaces, corresponding to the stationary mVN equation, is invariant under the full 
Lie sphere group). 

Remark. Diagonally-cyclidic surfaces have a natural projective counterpart, 
namely, the so-called isothermally-asymptotic surfaces for which the 3-web, formed 
by the asymptotic lines and Darboux's curves is hexagonal (Darboux's curves are 
the zero curves of Darboux's cubic form). This class of surfaces can be equivalently 
characterized by the existence of asymptotic coordinates where Darboux's cubic 
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form becomes isothermic. We refer to [|3lJ] for further discussion and references 
conserning isothermally-asymptotic surfaces (^-surfaces in the terminology of [|36|1 ). 
Isothermally-asymptotic surfaces are related to a different real reduction of the mVN 
equation. Here we present the details of it's derivation. Let M 2 be a surface in 
projective space parametrized by asymptotic coordinates x, y with the radius-vector 
r satisfying the equations 

^xx Q-Tx P^yi (a y\ 

r yy = qr x + bf y . 
With any surface (4.7) we associate the symmetric 2-form 

pq dxdy (4.8) 

and the conformal class of the Darboux cubic form 

pdx 3 + qdy 3 (4.9) 

which define "generic" surface M 2 uniquely up to projective equivalence and play a 
role similar to that of (1.1) and (1.2) in the Lie sphere geometry. Darboux's curves 
are the zero curves of the Darboux cubic form. The hexagonality conditions of the 
3-web formed by the asymptotic lines and the Darboux curves is equivalent to the 
existence of asymptotic parametrization x, y such that cubic form (4.9) becomes 
proportional to 

dx 3 + dy 3 , 

that is, to the condition p = q. An important subclass of isothermally-asymptotic 
surfaces are the proper affine spheres, for which the radius-vector r satisfies the 
equations 

xx p 1 x ~r jJi y 

/in Tl'T* Py iy 

1 yy — P 1 x p iy 

with p satisfying the Tzitzeica equation 

(lnp) xy = p 2 - i. 

In the case p = q the compatibility conditions of (4.7) reduce to 

(Px + ap+ \b 2 - b y ) x = \{p 2 ) y 

(Py + bp+la 2 -a x ) y = l(p 2 ) x (4.10) 

dy b X . 
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Equations (4.10) can be rewritten in the form 

by = Vx + ap + \b 2 - \W, W x = (p 2 ) y 

a x =p y + bp + \a 2 -|V, V y =(p 2 ) x (4.11) 

a y = f, b x = f. 
Crossdifferentiation of (4.11) gives the expressions for f X) f y 

fx = P yy + ap 2 + \pb 2 - \pW + bp y + af- (p 2 ) x 

(4.12) 

fy = Pxx + bp 2 + \pa 2 - \pV + ap x + bf - {p 2 ) y 
the compatibility conditions of which result in the stationary mVN equation 

PXXX - 3V Px - \pV X = Pyyy ~ 3W P y ~ \pWy 

W x = (p 2 ) y (4.13) 
Vy = (P 2 )x 

which can be reduced to that presented above by a complex change of variables 
p — > iU, W — > — W, V — > —V. Integrating the compatible systems (4.12), (4.11) 
and (4.7) for a given solution p, W, V of (4.13) we arrive at the explicit formula for 
the radius-vector r. 

Particular solutions of the stationary mVN equation (4.13) can be obtained by 
the ansatz 

w= 2 _Pyy_±(Py) 2 
3p 3\p) 

= 2p^ _ 1 fpA 2 
" 3 p 3\p) 
where p satisfies the Tzitzeica equation 



c 



(\np) xy = p H — , c = const. 

In this case the first equation is satisfied identically, while the last two become just 
two conservation laws of the Tzitzeica equation. If c ^ 0, then it can be normalized 
to — 1; the corresponding surfaces are the proper affine spheres. The case c = 
corresponds to the improper affine spheres whose affine normals are parallel and the 
radius- vector r satisfies the equations 

xx p I x ~ y> y 

ty — 'Y^y Py y 

yy P' x ply 

_^ ^ — * 

r xy = pi 
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where I is a constant vector (direction of the affine normal) and p satisfies the 
Liouville equation 

(\np) xy = p 2 . 

The fact that Tzitzeica's equation defines a subclass of solutions of the stationary VN 
equation was observed recently in ]37| (I would like to thank W. Schief for providing 
me with this reference). Our results give differential-geometric interpretation of this 
formal observation. 

Isothermally-asymptotic surfaces are known to possess Backlund transformations 
such that the initial and the transformed surfaces are two focal surfaces of a W- 
congruence, which preserves asymptotic lines and the Darboux curves. This can 
be the starting point for the modern approach to isothermally-asymptotic surfaces, 
their discretization in the spirit of [E35[], etc. 



5 Invariants of higher dimensional hypersurfaces 
in Lie sphere geometry 

In this section we announce several results on Lie sphere geometry on higher dimen- 
sional hypersurfaces, postponing the detailed proofs to a separate publication. 

Let M n be hypersurface with principal curvatures k % and principal covectors a/, 
so that the i-th principal direction of M n is defined by the equations uji = 0, j ^ 
i. It should be pointed out that generic hypersurface of dimension > 3 does not 
possess parametrization R l by the lines of curvature as in the 2-dimensional case. 
Differentiating covectors uj 1 and principal curvatures k % we arrive at the structure 
equations 

dJ = c) k u j A u k , (5.1) 

du = ky. 

Let also 

ds 2 = 9au l 
l 

be the first fundamental form of hypersurface M n . 

Theorem 2. A generic hypersurface M n (n > 3) is defined by quadratic form 

(14) 

h.ih.j 

— J — UJW 



'1 

(A* - ti) 2 



and the conformal class of cubic form (1.5) 

i 

uniquely up to Lie sphere equivalence. 
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As "generic" it is sufficient to understand a surface with k\ 7^ 0. This genericity 
assumption is essential since, for instance, there exist examples of Dupin hypersur- 
faces (that is, hypersurfaces with k\ = 0) which are not Lie-equivalent. Theorem 2 
is an analog of the corresponding theorem in projective differential geometry stating 
that hypersurface in projective space P n of dimension n > 4 is uniquely determined 
by the conformal classes of of its second fundamental form and Darboux's cubic 
form (see for the exact statements and further references). 

Remark 1. In case k\ 7^ cubic form (1.5) encodes all the information about 
the lines of curvature of hypersurface M n . Indeed, principal directions are uniquely 
defined as those directions where cubic form (1.5) reduces to the sum of pure 
cubes (without mixed terms). Moreover, principal directions are zero directions 
of quadratic form (1.4). However, this last condition does not define them uniquely 
as in the 2-dimensional situation. 

Remark 2. The invariant quadratic form (1.4) defines the invariant volume 
form, giving rise in the case n = 3 to the invariant functional 

uWu 3 , (5.2) 



(k 1 - k 2 )(k l - k 3 )(k 2 - k 3 ) 

the extremals of which should be called minimal hypersurfaces in Lie sphere ge- 
ometry in analogy with the 2-dimensional case. It does not look likely that this 
functional was investigated so far. 

Remark 3. In principle for n > 3 there exist additional Lie-sphere invariants 
besides those mentioned in Theorem 2, namely: 

1. The cross-ratios 

(k { - ti)(k n - k l ) 



(k n - - k l ) 

of any four principal curvatures. 
2. The covectors 



(5.3) 



1 ^ (5.4) 



{k l - - k l ) 

For instance, in case n = 3 we have three invariant covectors 

h l (h 2 — h 3 ) p/p_J,i\ 
n 1 — to 1 Vl 2 - 



(A; 1 - A; 2 )(A; 1 - A; 3 ) ' (k 2 - k l ){k 2 - k 3 ) 

U3(ul _ h.2\ 

q3 _ ft 3l fc fc ) 3 
(k 3 - k l )(k 3 - k 2 ) 

giving rise to the invariant quadratic form 
whose volume functional 

r r / 



15 



coincides with (5.2). 

3. Conformal class of the quadratic form 



(\ n-2 / \ n-2 

nC* 1 -*') u l2 + ...+g nn \ll(k n -k 1 )) oo n2 . (5.5) 

Lie-invariant class of hypersurfaces with conformally flat quadratic form (5.5) de- 
serves a special investigation. 

4. Differential dfl of the 1-form 



(5.6) 



In generic case k\ ^ objects (5.3)-(5.6) can be expressed through the forms 
(1.4) and (1.5). However they are important in the nongeneric situations, when 
some (or all) of k\ vanish so that (1.4) and (1.5) become zero. In particular, cross- 
ratios of principal curvatures play essential role in the study of Dupin hypersurfaces 
- see fl3~2"| , f33|. In this respect it seems interesting to understand the role of confor- 




mal class (5.5) and the 2-form dVt in the modern Lie-geometric approach to Dupin 
hypersurfaces. 

For hypersurfaces with nonholonomic nets of lines of curvature Theorem 2 leads 
to a nice geometric corollary, which we will discuss in the simplest nontrivial 3- 
dimensional case. Let us consider the structure equations (5.1) of the 3-dimensional 
hypersurface M 3 . Then there are only two possibilities: 

1. Holonomic case: all three coefficients c\\i cf 2 are equal to zero. Such 
hypersurfaces possess parametrization by the lines of curvature. 

2. Nonholonomic case: all three coefficients c^, c^, c 3 2 are nonzero. 

It immediately follows from the Peterson-Codazzi equations that for n = 3 inter- 
mediate cases are forbidden. In the nonholonomic case we can normalize covectors 
uj 1 , u 2 , u 3 in such a way that the structure equations assume the form 

doj 1 = au 1 A J 2 + bid 1 A u 3 + u 2 A lu 3 , 

duo 2 = puj 2 A uj 1 + quj 2 A uj 3 + u 3 A lu 1 , (5.7) 
du 3 = ru 3 A u 1 + sou 3 A u 2 + u 1 A u 2 . 



This normalization fixes uo l uniquely. As follows from the results of [24], in the 3- 
dimensional nonholonomic situation Peterson-Codazzi equations completely deter- 
mine quadratic form (1.4) and qubic form (1.5) through the coefficients a, b,p, q, r, s 
in the structure equations (5.7). Hence we can formulate the following result. 

Theorem 3. Nonholonomic 3-dimensional hypersurface M 3 is defined by its 
structure equations (5.7) uniquely up to Lie sphere equivalence. 

We can reformulate this result as follows: two 3-dimensional nonholonomic hy- 
persurfaces are Lie-equivalent if and only if there exists a point correspondence 
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between them, mapping the lines of curvature of one of them onto the lines of cur- 
vature of the other. Hence 3-dimensional nonholonomic hypersurface is uniquely 
determined by geometry of it's lines of curvature. 

This Theorem should remain valid for higher- dimensional hypersurfaces if we 
generalize the notion of " nonholonomicity" in a proper way (e.g. c* fc 7^ for all 
i 7^ j 7^ k which probably can be weakend). 



6 Reciprocal transformations of hydrodynamic type 
systems. Reciprocal invariants 

In this section we consider 2-component systems of hydrodynamic type 

u\ = vi(u)ui, i,j = 1,2 (6.1) 

which naturally arise in polytropic gas dynamics, chromatography, plasticity , etc. 
and describe wide variety of models of continuous media. The main advantage of the 
2-component case is the existence of the so-called Riemann invariants: coordinates, 
where equations (6.1) assume the diagonal form 

R] = X 1 (R)Rl, , s 

considerably simplifying their investigation. Any system (6.2) possesses infinitely 
many conservation laws 

h(R)dx + g{R)dt (6.3) 
with the densities h(R) and the fluxes g(R) governed by the equations 

d i g = \ i d i h, i = l,2 (6.4) 

(di = 8/ OR 1 ) which are completely equivalent to the condition h t = g x , manifesting 
closedness of 1-form (6.3). Crossdifferentiation of (6.4) results in the second-order 
equation 

dld<2h = \i dlh + \i (6.5) 

for the conserved densities of system (6.2). Thus conservation laws of system (6.2) 
depend on two arbitrary functions of one variable. Let us choose two particular 
conservation laws B(R)dx + A(R)dt and N(R)dx + M(R)dt and introduce new 
independent variables X, T by the formulae 

dX = Bdx + Adt , . 

dT = Ndx + Mdt ( > 
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which are correct since the right hand sides are closed. Changing from x,t to X, T 
in (6.2) we arrive at the transformed system 

R\ = A 1 (R)Rl , , 

R 2 = A 2 (R)R 2 X [bJ} 

where the new characteristic velocities A* are given by the formulae 

A«= X%B ~ A , i = l,2. (6.8) 

Remark. In principle one can apply transformation (6.6) directly to system (6.1) 
without rewriting it in Riemann invariants. In this case the transformed equations 
assume the form 

Ut = Vj{u)v? x , 

with the new matrix V given by 

V = (Bv - AE) (ME - Nv)~\ E = id. 

Transformations of type (6.6) are known as "reciprocal" and have been exten- 
sively investigated in [JIB] , |17| (see also |22| and |1| - pl| for further discussion). 
Following ]I8 |, [IS] we introduce the reciprocal invariants: 

the symmetric 2-form 



dlX ' d2X2 dR'dR 2 (6.9) 



and the differential 
of the 1-form 



(A 1 -A 2 ) 2 

dQ (6.10) 



f mx 2 d l X 1 \ . ( d&# d 2 X 2 \ 2 

n= {-d^ + xT^) dR + {-^xr + xJ^) dR (6 ' n) 

(Q itself is not reciprocally invariant). Note that both objects (6.9) and (6.10) do 
not change under the reparametrization of Riemann invariants R l — > ip 1 (R 1 ) ) R 2 — > 
<p 2 (R 2 ). 

Remark. In order to check the invariance of (6.9) and (6.10) under arbitrary 
reciprocal transformations it is sufficient to check their invariance under the following 
elementary ones: 

dX = Bdx + Adt, 
dT = dt, 

which changes only x and preserves t (under this transformation A* goes to A* = 
X l B - A) and 

dX = dt, 
dT = dx, 
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which transforms A* into A* = jj. The invariance of (6.9) and (6.10) under these 
elementary transformations can be checked by a direct calculation. Since any recip- 
rocal transformation is a composition of elementary ones, we arrive at the required 
invariance. 

It is quite remarkable that invariants (6.9) and (6.10) form a complete set in the 
following sense: if the invariants of one system can be mapped onto the invariants of 
the other one by the appropriate change of coordinates R\ than both these systems 
are reciprocally related and the corresponding reciprocal transformation (6.6) can 
be constructed explicitely (see |18[, |3| for the discussion). 



7 Hamiltonian systems 

System (6.1) is called Hamiltonian, if it can be represented in the form 

with the Hamiltonian operator and the Hamiltonian H = J h(u)dx. In this 

case the matrix v * is just the Hessian of the density h (for defmiteness we choose 
e l — 1), so that system (6.1) assumes the form 



u 1 \ ( /in hi2 \ ( u 1 
u 2 I { h 12 h 22 I { u 2 



(7.1) 



(here means a a ~ ). For systems (6.2) in Riemann invariants the necessary and 



sufficient condition for the existence of the Hamiltonian representation (7.1) is given 
by the following 

Lemma ||25|| . System (6.2) is Hamiltonian if and only if there exists flat diagonal 
metric ds 2 = gn^IfydR 1 + g 2 2(R)dR 2 such that 

d 2 l nv ^E= $^r, 

(7.2) 

«9il nv ^=^. 

Introducing the Lame coefficients Hi = H 2 = \[g 22 and the rotation coeffi- 

cients (3i2,(3 2 i by the formulae 

diH 2 = (3 12 H h d 2 Hi = /3 2 iH 2 , (7.3) 

we can rewrite the flatness condition of the metric ds 2 in a simple form 

diPw + d 2 (3 2 i = 0. (7.4) 

The coordinates m 1 ,^ 2 in (7.1) are just flat coordinates of the metric ds 2 , where it 
assumes the standard Euclidean form du l + du 2 . 
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Under reciprocal transformations (6.6) the metric coefficients gu transform ac- 
cording to the formulae 



(M - AW) 2 
(BM - ANf 



g « = 0htWf 7^, < = 1,2 (7.5) 



sec 



22] , [p3f] ), so that the transformed metric coefficients and the transformed 
characteristic velocities A* satisfy the same equations (7.2). It is important to em- 
phasize that reciprocal transformations do not preserve in general the flatness con- 
dition of the metric ds 2 and hence destroy the Hamiltonian structure. However, for 
any Hamiltonian system there always exist sufficiently many "canonical" reciprocal 
transformations preserving the flatness condition |f23 . 
Let us introduce the cubic form 

d^gndR 13 + d 2 \ 2 g 22 dR 2 \ (7.6) 

Using formulae (6.8) and (7.5) one can immediately check, that this qubic form 
is conformally invariant under reciprocal transformations: it aquires the multiple 
bm-an ' so ^ na ^ ^ ne zero curves °f (7-6) are reciprocally invariant. Hence with 
any Hamiltonian system we can associate besides the invariants (6.9) and (6.10) 
the reciprocally invariant 3-web of curves formed by coordinate lines R 1 = const, 
R 2 = const and the zero curves of qubic form (7.6) which are defined by the equation 



(d^gn^dR 1 + (d 2 \ 2 g 22 )^dR 2 = 0. 



A calculation similar to that in sect. 4 shows that invariant (6.10) is just the curvature 
form of this 3-web. 

Remark. It will be interesting to obtain explicit formulae for reciprocal invari- 
ants (6.9), (6.10) and (7.6) in the flat coordinates u % in terms of the Hamiltonian 
density h. 

As we already know the objects similar to (6.9), (6.10) and (7.6) arise in the Lie 



sphere geometry of surfaces. To clarify this point we recall the construction of [20 
[pl|j relating Hamiltonian systems and surfaces in the Euclidean space. 



8 Hamiltonian systems and surfaces in E 3 

Let us consider a 2-component Hamiltonian system (7.1) 




and apply the reciprocal transformation 

dX = Bdx + Adt, 
dT = dt, 
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where 



B 



u l2 + u 22 + l 



A = hiu 1 + h 2 u 2 - h 



(this is indeed an integral of system (7.1)). The transformed system assumes the 
form 



u l \ ( h n B - A h 12 B \ ( u l 



u 



hioB ho 2 B-A 



u 



To reveal geometric meaning of system (8.1) we introduce a surface M 2 in the 
Euclidean space E 3 (x 1 ,x 2 ,x 3 ) with the radius-vector 



( h-u 1 ^ \ 



r = 



x 

9 

X 



X 



ho 



u 2 — 

a B 



V 



A 
B 



(8.2) 



/ 



As one can verify by a straightforward calculation, the unit normal of the surface 
M 2 is given by 

( — \ 

1 B 



n 



B 



v i-i J 



Let us define the matrix w) by the formula 



dr Jiy , dn 

Geometrically is just the inverse of the Weingarten operator (shape operator) of 
the surface M 2 . Using the formulae for r and n we arrive at the following expression 
for the matrix Wj\ 

( huB-A h l2 B \ 
\ h l2 B h 22 B-A ) 

which coincides with that in (8.1). Hence the matrix of system (8.1) is just the 
inverse of the Weingarten operator of the associated surface M 2 . The characteristic 
velocities w % of system (8.1) are related to that of system (7.1) by the formula 



w i = XB - A, 



(8.3) 



and have geometric meaning of radii of principal curvatures of the surface M 2 . 
Moreover, the Riemann invariants of both systems (7.1) and (8.1) coincide and play 
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the role of parameters of the lines of curvature. Equations (8.1) can be equivalently 
represented in the conservative form 



n t = r x 



Some further properties of the correspondence (8.2) (in the general n-component 
case) were discussed in p0| , pl| , in particular: 



- commuting Hamiltonian systems correspond via formula (8.2) to surfaces with 
the same spherical image of the lines of curvature; 

- multi-Hamiltonian systems correspond to surfaces, possessing nontrivial de- 
formations preserving the Weingarten operator; 

- " canonical" reciprocal transformations, preserving the Hamiltonian structure, 
correspond to Lie sphere transformations of the associated surfaces; 

- the flat metric ds 2 defining Hamiltonian structure (see lemma in sect. 7) cor- 
responds to the third fundamental form of the associated surface. 

Since the correspondence between systems (7.1) and (8.1) is reciprocal, invariants 
(6.9), (6.10) and (7.6) coincide respectively with the symmetric 2-form 

dlwld2Uj2 dR 1 dR 2 , 



(w 1 — w 2 ) 2 

the skew-symmetric 2-form dQ, where 

n = ( d x d 2 w 2 + <W \ dRl + f d x d 2 w l + d 2 w 2 \ dR2 
\ d 2 w 2 w 1 — w 2 J \ diw 1 w 2 — w 1 J 

and the conformal class of the qubic form 

d^GndR 13 + d 2 w 2 G 22 dR 23 

where now R},R 2 are the parameters of lines of curvature, w 1 ^ 2 are the radii of 
principal curvatures and Gu, G 22 are the components of the third fundamental form 
of the associated surface M 2 . Since these objects preserve their form if we rewrite 
them in terms of principal curvatures k l and the components of the metric gu, they 
concide with the Lie sphere invariants of the surface M 2 . This provides remarkable 
differential-geometric interpretation of reciprocal invariants of hydrodynamic type 
systems. 

9 Reciprocal transformations and reciprocal in- 
variants of n-component systems 

Let us consider an n-component system of hydrodynamic type 

u] = vUu)u j x , i,j = l,...,n (9.1) 
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with the characteristic velocities A* and the corresponding left eigenvectors l l = (/*■) 
which satisfy the formulae 

k 

Introducing the 1-forms of = l^du^ (note that l l and oj 1 are defined up to rescailing 
ji _^ pty^ p 1 ^^ we can rewrite equations (9.1) in the equivalent exterior form 

of A (dx + X i dt) — 0, i = l,...,ra. (9.2) 

Differentiation of cj 1 and A* results in the "structure equations" of system (9.1): 

dJ = c) k u j A oj\ (9.3) 

dA* = \)uj j . (9.4) 

Systems in Riemann invariants are specified by the conditions Cj k = for any triple of 
indices % ^ j ^ k. Indeed, in this case the forms uj % satisfy the equations duj 1 Au l = 
for any i — 1, and hence can be normalized so as to become just ui l = dR\ In 
the coordinates R l equations (9.2) assume the familiar Riemann-invariant form 

Rl = \ i Rl i = l,..,n. (9.5) 

The exterior representation (9.2) is a natural analog of representation (9.5) which 
is applicable in the nondiagonalizable case as well. We emphasize that for n > 3 
Riemann invariants do not exist in general. 

Applying to (9.1) the reciprocal transformation 

dX = Bdx + Adt, 
dT = Ndx + Mdt, 

we arrive at the transformed equations 

U T — V j( u ) u X 

with the new matrix V given by 

V = (Bv - AE)(ME - Nv)-\ E = id 
or, in the exterior form, 

J A (dX + K l dT) = 

where 

A<= X ' B ~ A 



M -\ l N 

Hence the forms u l as well as the structure equations (9.3) do not change, while A* 
transform as in the 2-component case — see formula (6.8). 
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Remark. In the n-component case equations (6.4) for the densities and fluxes 
of conservation laws hdx + gdt assume the form 



X l h i} i — 1, n 



where Qi and hi are defined by the expansions 



dg = giiv\ dh = hiU 1 . 



In [081, 19 we introduced the following reciprocally invariant objects: 



1. The symmetric 2-form 



E TV^k? ^ ( 9 - 6 ) 



2. The skew-symmetric 2-form 
where 



dD, (9.7) 




x k A:: 



<-> = I '12 -TT-T7 1 I + • • • + l^ie^J ^ ( 9 - 8 ) 

(fl itself is not reciprocally invariant). Note that both objects (9.6) and (9.7) do 
not change if we reparametrize the 1-forms in the structure equations: uo l — > p l uj % . 
Objects (9.6) and (9.7) are the natural analogs of the corresponding invariants (6.9) 
and (6.10) in the 2-component case. However, for n > 3 the form Q depends only 
on the first derivatives of the characteristic velocities A' rather that on the second 
derivatives as in the 2-component case. 

In principle for n > 3 there exist additional reciprocal invariants, namely the 
1-forms 

(V-V)(A'*V <*^» 
as well as the cross-ratios 

(\ l - \i)(\ k - \ l ) 
(X k - X^>)(X i - X 1 ) 
of any four characteristic velocities. 



However, as follows from |]18 |, [19] , the structure equations (9.3) and the invari- 
ants (9.6), (9.7) in fact define generic system of hydrodynamic type uniquely up 
to reciprocal equivalence (under "generic" it is sufficient to understand genuinely 
nonlinear system, that is, a system with A* ^ for any i). 
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10 Appendix. Derivation of the mVN equation 
for diagonally- cyclidic surfaces 

Our aim is to show, that system (4.4) on p, k 1 , k 2 can be rewritten in terms of the 
Lie sphere density U defined as 



U 



, d^d^k 2 



(k 1 - k 2 y 



In view of the correspondence between surfaces in E 3 and Hamiltonian systems dis- 
cussed in sect. 8 the classification of diagonally-cyclidic surfaces can be equivalently 
reformulated as the classification of those Hamiltonian systems 

R\ = A 1 /?* 

E>2 _ \2 p2 

for which the flat metric ds 2 = gndR l2 + g22dR? 2 defining the Hamiltonian structure 
can be represented in the form 

e 2p e 2p 



911 = TTTT' 922 



d 1 \ 1 ' d 2 \ 2 ' 

Indeed, this is an immediate consequence of isotermicity of the cubic form. Equa- 
tions (7.2) and the flatness condition for the metric ds 2 imply the following system 
for p, A 1 , A 2 : 

Bin- B+ l 9id2>L. 
V\p — D -t- 2 Q ^ 2 , 



2 p - A + 2^TX^' 



10.1) 



^A 1 (d,B + |ft In + d 2 \ 2 [d 2 A + ^d 2 In |£) = (10.2) 

where A = jtzjt, B = (compare with (4.4)). Since 

U 2 



d 1 k 1 d 2 k 2 d 1 X 1 d 2 X 2 



(k 1 - k 2 ) 2 (A 1 - A 2 )' 

the desired equation for U will be obtained after we rewrite (10.1) - (10.2) in terms 

1 l„ d 2 x 2 

2 d^X 1 

dM = d x A - d 2 B, 



of U. Crossdifferentiating (10.1) and introducing k — \ In we obtain 



while (10.2) assumes the form 

d 2 {e k A) +d 1 {e~ k B) = 0. 
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Moreover, the equations for the characteristic velocities 

d 1 X 1 = Ue~ k (X 1 - X 2 ), <9 2 Ai = A(X 2 - A 1 ) 

d l X 2 = B(X l - A 2 ), d 2 X 2 = Ue^X 1 - A 2 ) 
give as the compatibility conditions two additional relations 

dxA = AB + U 2 -d 2 (Ue~ k ), 

d 2 B = AB + U 2 + d 1 (Ue k ). 

Thus the problem is reduced to that of excluding variables A,B,k and deriving the 
equation for U from the following much simpler looking system 

did 2 k = dxA - d 2 B, (10.3) 

d 2 (e k A)+d 1 (e- k B) = 0, (10.4) 

d ± A = AB + U 2 -d 2 (Ue~ k ), (10.5) 

d 2 B = AB + U 2 + d 1 (Ue k ). (10.6) 
To proceed further we introduce new variables m, n by the formulae 

-d 2 (Ue- k ) = Urn, d 1 (Ue k ) = Un, 

so that dik, d 2 k can be expressed as follows: 

d^k = + e~ k n, d 2 k = ^ + e k m. (10.7) 

The compatibility conditions of (10.7) with (10.3) (which now assumes the form 
d\d 2 k = U(m — n)) give the following equations for m, n: 



d 2 n = + e k {d 1 d 2 In U + U{m - n) + mn) 

dim = ^j-m — e~ k (did 2 laU — U(m — n) + mn). 
Rewriting (10.4) in the form 

e k d 2 A + e k A(^- + e k m) + e~ k d x B - e ~ k B(-^- + e~ k n) = 

and introducing the new variable F by the formulae 

e k d 2 A + e k A(^- + e k m) = F + -e k A 2 - -e~ k B 2 , 

e -k diB _ e -k B (_<W + e -* n ) = - F - \e k A 2 + \e~ k B 2 , 
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(10.8) 



we can express the derivatives of A and B as follows 

<9iA = AB + U 2 + Um, 

d 2 B = AB + U 2 + Un, 

d 2 A = -A{%f + e k m) + \A 2 - \e~ 2k B 2 + e~ k F, 

dxB = -B(&j?- - e- k n) - \e 2k A 2 + \B 2 - e k F. 
Let us introduce also the new functions G and H by the formulae 

d in = e 2k F + |e" fc (3n 2 + AnU + U 2 ) - 2n^§- - ZdtfJ + \e k H, 

d 2 m = e~ 2k F - \e k {3m 2 + AmU + U 2 ) - 2m^- - 3d 2 U + \e~ k G. 



(10.9) 



(10.10) 



The compatibility conditions of eqns. (10.9) give the following expressions for 
(h I '- d 2 F: 



d ± F = -F(2$- - e~ k n) + e 2k Um(U + 



m) 



+ f (G + 2e~ k F - e 2k (3m 2 + AmU + U 2 )), 

(10.11) 

d 2 F = -F{%f + e k m) + e~ 2k Un(U + n) 

-%(H + 2e k F + e~ 2k (3n 2 + AnU + U 2 )). 

The compatibility conditions of (10.8) and (10.10) give the expressions for d±G, d 2 H 
which can be represented as follows: 

d 1 (G + 2 9 -^-(°f) 2 ) = -3d 2 (U 2 ), 

(10.12) 

d 2 (H-2 d -^ + (^) 2 ) = 3d 1 (U 2 ). 

Finally, the compatibility conditions of (10.11) imply 

d 2 (GU 2 ) + d^HU 2 ) = 0, (10.13) 

so that equations (10.3)-(10.6) are reduced to (10.12)-(10.13). The desired stationary 
mVN equation results now after introducing V, W by the formulae 



All these calculations were performed with Mathematica and can be easily verified. 
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11 Concluding remarks 



Here we just list some of the unsolved problems. 

1. In our discussion of Lie-sphere invariants of surfaces (reciprocal invariants of 
hydrodynamic-type systems) the choice of coordinates R l plays a crucial role. In 
case of surfaces these are coordinates of the lines of curvature (Riemann invariants 
in case of hydrodynamic type systems). This choice is not accidental, since the lines 
of curvature are preserved by the Lie sphere group while Riemann invariants are 
preserved under reciprocal transformations. In fact only in these special coordinates 
do our invariants assume particularly symmetric and simple form. However, from 
the point of view of applications it is desirable to have a kind of invariant tensor 
formula, which will allow computation of these objects in an arbitrary coordinate 
system, for instance, in conformal parametrization in case of surfaces or in flat 
coordinates in case of Hamiltonian systems. 

2. In |]], [Hj Konopelchenko introduced dynamics of surfaces, governed by 



the mVN equation. In this approach the integrals of mVN correspond to certain 



functionals on surfaces, which were conjectured in ||12|| , | fL3|| , fl4| to be conformally 
invariant. In particular, the simplest quadratic integral of mVN corresponds to the 
Willmore functional. This conjecture was proved recently in ||15|| . Since functional 
(1.3) is conformally invariant (indeed, it is invariant under the full group of Lie 
sphere transformations which contains conformal group), it would be interesting to 
understand its relationship to the mVN hierarchy. 

3. It seems to be an interesting problem to describe the class of surfaces, for 
which evolutions (3.5), (3.6) preserve the lines of curvature. 
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